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An up to date review of Standard Big Bang Nucleosynthesis predictions vs the astrophysical estimates of light 
nuclei abundances is here presented. In particular the analysis reports the expected ranges for baryon fraction 
and effective number of neutrinos as obtained by BBN only. 



1. Introduction 

The nucleosynthesis taking place in the primor- 
dial plasma plays a twofold role: it is certainly one 
of the observational pillars of the hot Big Bang 
model, being indeed known simply as "Big Bang 
Nucleosynthesis" (BBN); at the same time, it pro- 
vides one of the earliest direct cosmological probe 
nowadays available, constraining the properties of 
the universe when it was a few seconds old, or 
equivalently at the MeV temperature scale. 

The basic framework of the BBN emerged in 
the decade between the seminal Alphcr-Bcthe- 
Gamow (known as afij) paper in 1948 [1] and the 
essential settlement of the paradigm of the stellar 
nucleosynthesis of elements heavier than 7 Li with 
the B 2 FH paper [2]. This pioneering period — an 
account of which can be found in [3] — established 
the basic picture that sees the four light-elements 
2 H, 3 He, 4 He and 7 Li as products of the early fire- 
ball, and virtually all the rest produced in stars 
or as a consequence of stellar explosions. 

In the following decades, the emphasis on the 
role played by BBN has evolved significantly. In 
the simplest scenario, the only free parameters 
in primordial nucleosynthesis are the baryon to 
photon ratio r\ (equivalently, the baryon density 
of the universe) and the neutrino chemical po- 
tentials, /j„ q . However, only neutrino chemical 
potentials larger than r\ by many orders of mag- 



nitude have appreciable effects. This is why the 
simple case where all fx Va 's are assumed to be 
negligibly small (e.g., of the same order of rf) 
is typically denoted as Standard BBN (SBBN). 
However, it is usual to add to 77, as an additional 
parameter, the so-called effective number of neu- 
trinos, Acff, which measure the amount of rela- 
tivistic d.o.f. at the time of BBN. 

2. Standard BBN theoretical predictions 
versus data 

The goal of a theoretical analysis of BBN is to 
obtain a reliable estimate of the model param- 
eters, once the experimental data on primordial 
abundances are known. In this Section we will 
consider only the case of standard BBN, where 
the only two free parameters are the value of 
the baryon energy density parameter f^ft 2 (or 
equivalently the baryon to photon number den- 
sity, 77) and possibly, a non standard value for the 
relativistic energy content during BBN. The lat- 
ter, after annihilation can be parameterized in 
terms of the effective number of neutrinos 
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As it is well known, iV c ff can differ from zero 
also via the non thermal equilibrium terms which 
characterize the relic neutrino distributions. We 
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Tabic 1 

N e g and AY P obtained for different cases, with 
and without neutrino oscillations, as reported in 
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Figure 1 . Frozen distortions of the flavor neutrino 
spectra as a function of the comoving momentum, 
for the best fit solar and atmospheric mixing pa- 
rameters. R is the scale factor. In the case where 
we allow for $13 7^ consistently with present 
bounds (blue dotted lines), one can distinguish 
the distortions for and v T (middle and lower, 
respectively). From 0]. 
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Table 2 

Comparison of the exact BBN results with a 
fixed- AN c s approximation. From [4]. 



show in Fig. [T] the asymptotic values of the fla- 
vor neutrino distribution, both without oscilla- 
tions and with non-zero mixing. The dependence 
of the non-thermal distortions in momentum is 
well visible, which reflects the fact that more en- 
ergetic neutrinos were interacting with e^ 1 for a 
longer period. Moreover, the effect of neutrino 
oscillations is evident, reducing the difference be- 
tween the flavor neutrino distortions. Fitting for- 
mulae for these distributions are available in [?]. 
In Table [T] we report the effect of non instanta- 
neous neutrino decoupling on the radiation en- 
ergy density, N e g, and on the 4 He mass fraction. 
By taking also into account neutrino oscillations, 
one finds a global change of AY p ~ 2.1 xlO -4 
which agrees with the results in [5] due to the 
inclusion of QED effects. Nevertheless the net ef- 
fect due to oscillations is about a factor 3 smaller 
than what previously estimated, due to the fail- 
ure of the momentum-averaged approximation to 
reproduce the true distortions. 



In particular, the N e g reported in Table[T]is the 
contribution of neutrinos to the whole radiation 
energy budget, but only at the very end of neu- 
trino decoupling. Hence, not all the AN e ff there 
reported will be really contributing to BBN pro- 
cesses. In order to clarify this subtle point, we 
report in Table [2] the effect on all light nuclides, 
of the non instantaneous neutrino decoupling in 
the simple scenario of no neutrino oscillation, and 
compare this column with the simple prescrip- 
tion of adding a fix AN c g = 0.013 contribution 
to radiation. Even though Y p is reproduced (by 
construction), this is not the case for the other 
nuclear yields. Similar analysis have been re- 
cently presented by various groups, which might 
be slightly different depending on the adopted 
values of Y p and/or 2 H/H experimental determi- 
nation, see e.g. [6171819110111112113114115116] . 

In the minimal BBN scenario the parame- 
ters reduces to the baryon density only, since 
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Figure 2. Values of the primordial abundances 
as a function of rjio, calculated for AN c g = 0. 
The hatched blue bands represent the experimen- 
tal determination with 1 — a statistical errors 
on Y p , 2 H, and 7 Li, while the red band is the 
upper bound obtained in Rcf. [T7]. Note that 
for high value of ryio all 7 Li comes from 7 Bc ra- 
dioactive decay via electron capture. The verti- 
cal green band represents WMAP 5-ycar result 
n B h 2 = 0.02273 ± 0.00062 [18]. 



A-ZVoff is just produced by the above non ther- 
mal distortions. Fig. [2] shows the dependence on 
?7io = r\ ■ 10 10 of the final value of the primordial 
yields, calculated using the nucleosynthesis code 
PArthENoPE [19] , along with the experimental val- 
ues of the abundances and their corresponding 
uncertainties, as discussed in details in Rcf. [20]. 
Just to summarize, according to the analysis of 
[20] the inferred primordial abundances result to 
be 



2 H/H = (2.87t°; 22 ) x 10~ 5 
3 He/H< (1.1 ±0.2) x 10~ 5 



Y p = 0.247 ± 0.002 stat ± 0.004 syst 
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To get confidence intervals for 77, one construct a 
likelihood function 

£(77)ocexp(-x 2 (r/)/2) , (6) 
with 

x 2 ( v ) = Y^[Mv)-x° bs }w i MlxM-x? s } -(7) 

ij 

The proportionality constant can be obtained by 
requiring normalization to unity, and Wij(r)) de- 
notes the inverse covariance matrix, 



Wijirj) = [a 2 tj + o? iexp 6ij + cr 2 Jother ] 
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where <7y and a 



i.exp represent the nuclear rate un- 
certainties and experimental uncertainties of nu- 
clide abundance Xj, respectively (we use the nu- 
clear rate uncertainties as in Ref. [H]), while by 



a 



ij. other 



we denote the propagated squared error 
matrix due to all other input parameter uncer- 
tainties (r n , Gn, etc.). 

We first consider 2 H abundance alone, to illus- 
trate the role of deuterium as an excellent bary- 
omcter. In this case the best fit values found 
are tt B h 2 = 0.021 ± 0.001 (r) 10 = 5.7 ± 0.3) at 
68% C.L. , and n B h 2 = 0.021 ± 0.002 at 95% 
C.L. . A similar analysis can be performed us- 
ing 4 He. In case only statistical error is consid- 
ered we get n B h 2 = 0.021™^ ( mo = 5.7±\i) 



0.02lt°;°J° at 95% 



at 68% C.L. , and Q B h 2 
C.L. . Fig. [3] shows the two likelihood pro- 
files, which nicely agree. When accounting for 
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Table 3 

The theoretical values of the nuclear abundances 
for some value of VL B h 2 . 
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Figure 3. Likelihood functions for 2 H/H (narrow) 
and Y p (broad) , using only the statistical error for 
the 4 Hc measurement. 



the largest possible systematic error on Y p , the 
determination of f2s/i 2 becomes even more dom- 
inated by deuterium. In any case, the result 
is compatible at 2-a with WMAP 5-ycar result 
n B h 2 = 0.02273 ± 0.00062 [TH]. The slight dis- 
agreement might have some impact on the de- 
termination from CMB anisotropics of the pri- 
mordial scalar perturbation spectral index 
noticed in [55], where the BBN determination of 
fish 2 from deuterium is used as a prior in the 
analysis of the five year data of WMAP. 

In Table[3]we report the values of some relevant 
abundances for three different baryon densities, 
evaluated using PArthENoPE [19] . Notice the very 
low prediction for 6 Li and that, for these values 
of baryon density, almost all 7 Li is produced by 
7 Be via its eventual electron capture process. 

If one relaxes the hypothesis of a standard num- 
ber of relativistic degrees of freedom, it is possi- 
ble to obtain bounds on the largest (or smallest) 
amount of radiation present at the BBN epoch, 
in the form of decoupled relativistic particles, or 



Figure 4. Contours at 68 and 95 % C.L. of the 
total likelihood function for deuterium and 4 He 
in the plane (ris/i 2 ,7V ff). The bands show the 
95% C.L. regions from deuterium (almost verti- 
cal) and Hclium-4 (horizontal), neglecting possi- 
ble systematic uncertainty on Y p . We also show 
the 95 % C.L. allowed region from Y p with sys- 
tematic error included (dashed lines). 



non standard features of active neutrinos. Fig. 
2] displays the contour plot of the total likeli- 
hood function, in the plane (Ob/i 2 ,^^), beau- 
tifully centered extremely close to the standard 
value N e s = 3.0. After marginalization one gets 
n B h 2 = 0.021 ± 0.001 and N cff = 2.97 ± 0.14 
at 68% C.L., and Q B h 2 = 0.021 ± 0.002 and 

N cS = 2.971EJ' 2 ? at 95% CL - ■ Notc that in 
this case we are using for Y p the statistical un- 
certainty only. More conservatively, if one also 
considers the systematic uncertainty on Y p , the 
allowed range for N c g becomes broader, 7V c ff = 
3.0±0.3 stat (2a)±0. 

"Jsyst ■ 
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3. Conclusions 

The "classical parameter" constrained by BBN 
is the baryon to photon ratio, 77, or equivalcntly 
the baryon abundance, fish 2 . At present, the 
constraint is dominated by the deuterium deter- 
mination, and we find Q, B h 2 = 0.021 ± 0.001(1 
a). This determination is consistent with the up- 
per limit on primordial 3 Hc/H (which provides 
a lower limit to rf), as well as with the range 
selected by He determinations, which however 
ever neglecting systematic errors provides a con- 
straint 5 times weaker. The agreement within 
2 a with the WMAP determination, Q B h 2 = 

0. 02273 ± 0.00062, represents a remarkable suc- 
cess of the Standard Cosmological Model. On 
the other hand, using this value as an input, a 
factor ~ 3 discrepancy remains with 7 Li determi- 
nations, which can hardly be reconciled even ac- 
counting for a conservative error budget in both 
observations and nuclear inputs. Even more puz- 
zling are some detections of traces of 6 Li at a level 
far above the one expected from Standard BBN. 
Both nuclides indicate that either their present 
observations do not reflect their primordial val- 
ues, and should thus be discarded for cosmo- 
logical purposes, or that the early cosmology is 
more complicated and exciting than the Stan- 
dard BBN lore. Neither a non-standard num- 
ber of massless degrees of freedom in the plasma 
(parameterized via N e g) or a lepton asymmetry 
£ e (all asymmetries assumed equal) can recon- 
cile the discrepancy. Current bounds on both 
quantities come basically from the He measure- 
ment, N e ft = 3.0 ± 0.3 sta t (2 a) ± 0.3 syst and 
£ e = 0.004 ± 0.017*** (2 a) ± 0.017 syst . 
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